Hole dynamics and photoemission in a t-J model for SrCu2(B03)2 
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The motion of a single hole in a t — J model for the two-dimensional spin-gap compound 
SrCu2 (B03)2 is investigated. The undoped Heisenberg model for this system has an exact dimer 
eigenstate and shows a phase transition between a dimerized and a Neel phase at a certain ratio of 
the magnetic couplings. We calculate the photoemission spectrum in the disordered phase using a 
generalized spin-polaron picture. By varying the inter-dimer hopping parameters we find a cross- 
over between a narrow quasiparticle band regime known from other strongly correlated systems and 
free-fermion behavior. The hole motion in the Neel-ordered phase is also briefly considered. 



Since the discovery of high-temperature superconduc- 
tivity, doped antiferromagnets (AF) have been stud- 
ied intensively. The pseudo-gap behavior observed in 
the high-T c cuprates has stimulated great interest in 
systems with spin gaps. Several new one- and two- 
dimensional spin gap systems have been found ex- 
perimentally. These materials are characterized by a 
disordered singlet ground state and a finite gap to 
all spin excitations. Some of the compounds which 
have two-dimensional (2d) character include the coup! 
spin ladder system^ SrCu 2 3 lil, CaV 2 O 5 0, (V0 2 )P 2 7 l 
Cu 2 (C5Hi 2 N 2 ) 2 CLia, and the plaquette RVB system, 
CaV 4 Ogla. 

Recently the two-dimensional spin gap system 
SrCu 2 (B03) 2 has been found by Kageyama et alE. It 
has a spin-singlet ground state with a spin gap ~ 30 K. 
The substance has additional interesting features, e.g., 
the high-field magnetization was observed to have two 
plateaus at 1/4 and 1/8 of the full moment. Recent 
workQ suggests that the underlying physics can be under- 
stood on the basis of a two-dimensional S = h Heisen- 
berg model with antiferromagnetic nearest-neighbor (Ji, 
on links A) and next-nearest-neighbor (J 2 , on links B) 
couplings on the lattice shown in Fig. pi. The nearest- 
neighbor bonds A define a unique singlet covering of the 
lattice. The Heisenberg model corresponding to Fig. |l| 
(with Ji, J 2 ) is in fact topologically equivalent to the 
model considered by Shastry and Sutherland!! For this 
model the singlet product state forms an exact eigen- 
state of the Hamiltonian at all couplings, and is the 
ground-state in a region where the nearest-neighbor cou- 
pling J\ dominates. On the other hand, for J\ — > the 
system becomes equivalent to the 2d square lattice AF 
(with nearest-neighbor coupling J 2 ) which has a Neel- 
ordered ground state. The Shastry-Sutherland model can 
be complemented by a coupling J3 on the links C; the 
singlet product state is an eigenstate of this generalized 
model, too. For J 2 = J3 the total spin on each A bond 
is conserved, i.e., there is a macroscopic number of con- 
served quantities, and each eigenstate is characterized by 
the number and positions of triplets. 

The Heisenberg model of Fig. 1 has been studied in 
Rcfs. [7]|]-[ll] using exact diagonalization, Schwinger bo- 
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FIG. 1. Lattice structure of the Cu spins in SrCu2(BOs)2, 
with the three different exchange couplings and the crystallo- 
graphic axes a,b. The dash-dotted lines denote a unit cell. 



son and series expansion methods. For J 2 = J3 it can 
be mapped onto a spin-1 model which shows a first order 
transition between the singlet state and a Neel-ordered 
state at J2/J1 = 0.4296. For J 2 7^ J 3 the situation is 
less clear. For J3 = the numerical results provide ev- 
idence for a weak first order singlet-Necl transition at 
J2/J1 = 0.70 ± 0.01. The disordered SrCu 2 (B0 3 ) 2 com- 
pound lies probably close to the transition line to a Neel 
state, which explains the unusual temperature depen- 
dence of magnetic properties. A fit of the susceptibility 
obtained from the-.model (with J± and J 2 only) to the 
experimental dataQ leads to the estimates of J\ = 100 
K and J2/J1 = 0.68. However, the actual ratio of the 
couplings J3/J2 is not known. 

To include the doping degree of freedom into the 
Shastry-Sutherland model we consider a standard t — J 
model on the SrCu 2 (B03) 2 lattice: 

n = E (-ti(&l6 jtr + h.c.) + J 1 Si-S j ) 
+ E (~*a (4A CT + h.c.) + J 2 St ■ S k 

(i,k)eB 

+ E ("*3 (4Ao + h.c) + J 3 S, • S,) . (1) 

(i,i)ec 
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The electron operators c ia exclude double occupancies. 
We have included hopping ts and interaction J3 along the 
C bonds. If the t — J model is derived as strong-coupling 
limit of a Hubbard model on the SrCu2(BOa)2 lattice 
with on-site repulsion U, the ratio of the parameters is 
given by 

t 2 t 2 t 2 11 

_1 — _ _3_ _ _ (0) 

Jx ~ J2 ~ J3 ~ 4 ' u 

Although the compound SrCu 2 (B0 3 )2 has to our 
knowledge not been doped so far, the study of the hole 
dynamics in this environment is an interesting and chal- 
lenging question. Furthermore it may be possible that 
finite doping leads to the formation of hole pairs and 
eventually to superconductivity. In this paper we shall 
discuss the dynamics of a single hole in an otherwise half- 
filled system using a generalized spin-polaron picture. 
The one-hole spectral function for this case corresponds 
directly to the result of an angle-resolved photoemission 
experiment on the undoped compound which may give 
important information on the electronic correlations and 
the exchange constants in SrCu2 (603)2- We will show 
that the ratio t% /t% tunes a cross-over between a narrow- 
band quasiparticlc behavior and a regime where a free- 
fcrmion peak dominates the spectrum. The narrow-band 
behavior found here is similar to other 2d AF systemsEj, 
it originates^am the motion of a hole dressed with spin 
fluctuationsEra. 

To investigate the hole motion we consider a one- 
particle Green's function describing the creation of a sin- 
gle hole with momentum k at zero temperature: 

G(k, W ) = £«|cL-^c kff |<) (3) 

where z is the complex frequency variable, z = lo + irj, 
77 — > 0. IV'oO i s the ground state of undoped system, i.e., 
with N electrons on N lattice sites. 

To describe the dimerized phase we employ a bond 
operator representatio n! 1 ^I 17 ! for the spins on the nearest- 
neighbor bonds. For each A bond containing two S = 1/2 
spins we introduce bosonic operators for creation of a 
singlet and three triplet states out of the vacuum |0): 

* f |0) = ^(1 UH it)), 410) = ^(| TTH U», 41°) = 

75(1 TT) + I ID), 410) - ^(1 U) + I IT)), where the 

constraint s^s + J2 a Ofa = 1 nas *° ^ e i m P ose d on each 
bond to restrict the possible states to the physical Hilbert 
space. The original spins are related to the new boson 

basis operators by S" 2 = §(i s ^ Q 4 s — itap-ytpt-y) ■ 
The Hamiltonian of the Shastry-Sutherland model writ- 
ten in terms of the bond operators contains no terms 
which create triplet excitations from a state containing 
only singlets which means that the singlet product state 
\(j>o) = rii s !|0) is an exact eigenstate of the undoped 
system at all couplings. 

If we remove one electron from an A bond, a single- 
hole state on this bond is created. We introduce fermionic 



operators for bonding (symmetric) and antibonding (an- 
tisymmetric) states of one electron (or hole) on an A 
bond: 

at iCT |0) = i(c^ + 4 iCT )|0), 

a^|0) = i=(c t liCT -cL)|0). (4) 

The operators c\ 2 a create an electron with spin a on 
one of the two sites of an A bond. Hopping via t\ leads 
to an 'on-bond' energy of ±ii for a a and a g , respectively. 
The other hopping terms in TL permit the hole to hop 
between bonds; the exchange interaction gives rise to spin 
fluctuations on neighboring bonds. The full interaction 
Hamiltonian in terms of s, t, and a operators can be 
easily found by calculating all possible one-hole matrix 
elements of the initial Hamiltonian TL (Q) (see e.g. Ref. 

The evaluation of the Green's function (H) is done us- 
ing the Mori-Zwanzig projection technique. The set of 
dynamic variables is constructed from generalized path 
operators tM3 which here create strings of triplet exci- 
tations attached to the hole. Details of the calculational 
procedure can be found in Refs. In the disordered 

phase of the Shastry-Sutherland model the evaluation of 
the matrix elements is simplified by the fact that the 
undoped ground state does not contain background spin 
fluctuations, i.e., no cumulant expectation values are in- 
volved (cf. Ref. gfj). In the present calculations we have 
employed up to 1800 dynamic variables with a maximum 
path length of 3. The neglect of the self-energy terms 
leads to a discrete set of poles for the Green's functions, 
so the present approach cannot account for lincwidths. 
In all figures we have introduced an artificial linewidth 
of 0.2fi to plot the spectra. 

Before turning to the hole dynamics it is worth men- 
tioning that even a static vacancy (equivalent to the limit 
t — > 0) has non-trivial consequences: The pure singlet 
state with one spin removed, Ci a \<Po)j is 110 longer an 
eigenstate of TL. The unpaired spin leads to triplet exci- 
tations in its neighborhood ('screening cloud'), these are 
spatially confined as long as the spin gap is non-zero. 

Now we consider the case of non-zero hopping. We 
start with the 'symmetric' choice t 2 — £3, J2 — J3, which 
allows to obtain several results analytically. J2 = J?, 
implies that triplets are strictly localized in the absence 
of the hole, so any triplet created by hole hopping re- 
mains on its bond until the hole returns. t 2 = £3 re- 
duces the number of possible inter-bond hopping pro- 
cesses. The only non-zero hopping matrix elements are 
\(a a s\H\sa a )\ = \(a s t a \H\t a a s )\ = \(a s s\H\t a a a )\ = t 2 
where \XY) = xJyJ\Q) is a short-hand notation for 
a state of two neighboring bonds i and j. It fol- 
lows that a hole in the antisymmetric state a a can 
freely propagate in a singlet background without emis- 
sion of triplet fluctuations (direct hopping), i.e., |^ a .k) = 
J2i exp(ikRi)aj l jS;|0o) is an exact eigenstate of TL. In 
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contrast, a hole being in the symmetric state a s always 
creates a triplet (and converts into a a ) when it hops to a 
neighboring singlet bond. This in turns means that any 
one-hole eigenstate of Ti. which contains components with 
a s also involves triplet excitations. Since the triplets are 
localized and can only be created/removed by the hole, 
the state itself is localized, and contributions from such 
eigenstates to the spectrum are momentum-independent 
(non-dispersive) . 



a) t } /t 2 = 1 



b) t 3 /t 2 = 1/2 




c) I 3 /! 2 = 1/4 



d) r 3 = 




— (Jt/2,Jt/2) 
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FIG. 2. Evolution of the one-hole spectral function 
— Im G(k, uS) under a change of the hopping ratio %/t2- a) 
'Symmetric' case t^/t2 = 1, J2/J1 = 0.4, b) tz/t?. = 1/2, 
J2/J1 = 0.5, c) ta/ts = 1/4, J2/J1 = 0.6 d) ta = 0, 
J2/J1 = 0.68. The other parameters are chosen according 
to eq. (Q) with U/ti = 4. The energies are measured in units 
of ti relative to the energy of a localized hole. The ratio of the 
magnetic couplings is chosen tensJace the system close to the 
boundary of the singlet phaseu'til, i.e., the spin gap in units 
of Ji is nearly equal in the four cases. 

The calculated spectrum (Fig. ||a) is therefore eas- 
ily understood: It shows two dispersing bands which 
correspond to hole hopping in the antisymmetric state 
a Q with effective hopping amplitude ti — £3 through 
the square lattice of rungs (two bands arise from the 
fact that the unit cell contains two rungs). All other 
contributions are localized and involve symmetric hole 
states. The corresponding wave functions can be mod- 
eled by a particle moving in an attractive potential cen- 



tered at a single site i. A reasonable approximation 
for the states contributing to the spectrum is (uaj { Si + 

Y,n, a v R at a,i+ii s i+R t i a s i)\ < l > o) where u,v are some coeffi- 
cients (typically \u\ ^> |wr| and wr rapidly decaying with 
distance). Components with more than one triplet have 
coefficients much smaller than the one-triplet coefficients 
wr. Variation of the model parameters (keeping t% = £3, 
J2 = J3) only results in small changes in the spectrum 
of Fig. |2|a since the dominating bands are determined 
by t% only. Note, however, that varying t\/t2 shifts the 
dispersing bands w.r.t. to the localized peaks, i.e., it can 
induce a level crossing at momentum (0,0), so t\ » t 2 
eventually leads to a localized one-hole ground state. 

Having understood the special situation at t 2 — £3, 
Ji = 'hi w e turn to the general case. We assume J3 < J 2 
[and £3 < t% because of (0) ] since the model behavior is 
symmetric with respect to the interchange of the B and 
C bonds. In the following we discuss the t and J pa- 
rameters separately, but we keep in mind that they are 
usually connected by the relation (g). For J 2 7^ J3 the 
triplets are no longer completely localized, however, as is 
known from the undoped model, triplet hopping does not 
occur up to 6th order of perturbation theory in J 2 , and 
so the triplet dispersion is weak. Consequently, J 2 7^ J3 
introduces a weak dispersion into the 'background fea- 
tures' of the spectrum. The effect of t 2 7^ £3 on the 
spectrum is more pronounced: Hopping of the antisym- 
metric hole state now also emits triplets, i.e., IVVk) is 
no longer an exact eigenstate of Ti. Therefore the well- 
defined dispersing bands of Fig. ||a mix with the back- 
ground features. Spectral weight is transferred to the 
bottom of the spectrum which is easily understood from 
the fact that any triplet excitation increases the energy 
by at least the gap size. Figs. ||a - ||d show the evo- 
lution of the spectrum when going from ts/t<z = 1 to 
£3/^2 = 0; we have fixed U/ti = 4 and chosen the re- 
maining parameters such that the size of the spin gap is 
(approximately) preserved. For <3 = J3 = 0, Fig. ||d, 
we arrive at a situation with a narrow band at the bot- 
tom of the spectrum with dispersion minimum at (0,0); 
this band corresponds to the motion of a hole surrounded 
by triplet fluctuations (spin polaron) . At higher energies 
a weak background is visible, it arises from excited po- 
laron states. In addition, we note that smaller relative 
hopping strengths t/J suppress triplet fluctuations and 
reduce the high energy background, whereas larger val- 
ues of t/J lead to an incoherent spectrum since the size 
of the spin polaron increases. 

To make contact with possible experiments we exam- 
ine briefly the properties of the pronounced bands at the 
lower edge of the spectrum which would be visible in a 
photoemission spectrum. The experimentally measured 
bandwidth can be used to determine the ratio t/J, or 
equivalently, the on-site repulsion U. For the (experi- 
mentally unrealistic) case of t% = £3, J2 — J3, two bands 
should be observed as in Fig. ||a. The width of each of 
the bands is given by 4i 2 since the dispersion is the one 
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FIG. 3. Comparison of hole QP bandwidths in units 
of the (inter-dimer) in-plane hopping amplitude. Solid: 
SrCu 2 (B0 3 ) 2 lattice, t 3 /*i = 0.82, J 2 /Ji = 0.68, t 3 = J 3 = 0, 
present calculation. Dashed: Square lattice, obtained by the 
methods used in this paper. Circles: Square lattice SCBA 
results from Ref. 14l 



of a free fermion. More likely, the material has t% <C t-z, 
■h -C J2, which corresponds to Fig. |c ord. As a guide 
we plot in Fig. || the bandwidth for £3 = J3 = as 
function of J\/t\. In contrary to the well-known square- 
lattice case (shown for comparison) the bandwidth is fi- 
nite for J — > 0. The reason is the possibility of direct 
hopping (without emission of spin excitations); this is 
not the case in a square lattice antiferromagnetic back- 
ground where hopping always creates spin defects which 
have to be removed by exchange processes. For tj J — > 
the bandwidth saturates (in units of t^) at a non-trivial 
value which arises from the finite cloud of triplet excita- 
tions around a static hole. (Of course, for J2 = J3 = 
this effect is absent, and the bandwidth becomes 2t2 in 
the limit of t/J ^ 0.) 

We have also calculated the spin correlations near the 
mobile hole. Except for the case corresponding to Fig. 
||a which shows no inter-dimer correlations in the ground 
state, the hole always introduces antiferromagnetic cor- 
relations between different dimers in its vicinity. How- 
ever, for the possibly relevant values (J2/J1 = 0.68, 
U/ti ~ 2 — 10) the dimerization is strong and the po- 
laron size can be estimated to be smaller than 3 lattice 
spacings. 

We briefly mention that the present analysis can be 
extended to the Neel-ordered phase of the Shastry- 
Suthcrland model as well. Therefore a condensate of one 
type of triplets is introduced by|-a. proper transformation 
of the basis states on each bonded; the condensation am- 
plitude can be extracted from series expansion resultslij. 
The ground state of the undoped system is obtained by 
an expansion around a modified product state (which 
is the Neel state in the case of dominating J2); here of 
cause the ground state contains fluctuations around the 
product state. The spin polaron consists of spin devia- 
tions from the undoped background state in the vicinity 
of the hole. The results show that with increasing anti- 
ferromagnetic correlations the band minimum is shifted 
from (0,0) to (71", tt) (this is equivalent to (tt/2,it/2) in 
the Brillouin zone of the square lattice defined by the B 



bonds); the bands approach the shape known from the 
square-lattice antiferromagnet. 

Summarizing, in this paper we have studied the one- 
hole dynamics in a t — J model for the 2d spin gap mate- 
rial SrCu2(B03)2- Using a generalized spin-polaron con- 
cept together with an expansion around a singlet product 
state we have calculated the one-hole spectral function. 
It shows an interesting cross-over from free-fermion be- 
havior to correlated behavior under a variation of the 
ratio ta/i3 of 2nd and 3rd- nearest-neighbor coupling. At 
£3 = ^2 the main contributions to the spectrum can be 
described by the hopping of the bare hole between the 
singlet rungs with a tight-binding dispersion. In con- 
trast, for £3 <§C <2 one finds a narrow band at the bottom 
of the spectrum and a background at higher energies; 
this structure can be attributed to the motion of a hole 
dressed with spin fluctuations (similar to other strongly 
correlated systems). We add that very large t/J leads 
to one-hole ground states with higher spin and ferromag- 
netic correlations in spirit of the Nagaoka effect, this has 
not been considered here. The hole dynamics at finite 
doping and the possibility of hole pairing are interesting 
subjects of future research. 

Financial support from the DFG (VO 794/1-1) is 
gratefully acknowledged. 
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